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Non-invasive  Phosphorus -31  MR  spectral  characterization  of 
in  vivo  breast  tissue  anomalies  using  pattern  recognition 
and  artificial  intelligence  (AI)  techniques 


I .  Introduction : 

It  is  highly  desirable  to  develop  a  non-invasive 
spectroscopic  and  pattern  recognition  technique  that  can 
detect  and  reliably  interpret  images  or  spectral  data  from 
small  volumes  of  the  breast.  Due  to  the  pervasive  nature 
of  breast  cancer  in  society  today,  and  the  consequent  need 
of  a  highly  accurate,  early  diagnostic  tool,  this  is  a 
very  timely  proposal  that  could  have  a  significant  impact 
on  women's  health.  Patient  Rotating  Delivery  of 

Excitation  Off-resonance  (RODEO)  MRI  data  has  been 
obtained  from  Dr.  Diana  Lindquist  at  the  University  of 
Arkansas  for  Medical  Sciences.  These  patients,  which 
flagged  suspicious  regions  in  breast  tissue,  have 
undergone  needle  biopsies  from  these  suspect  regions  for 
pathological  examination.  With  the  patient's  permission. 
Dr.  Lindquist  obtained  P-31  MR  scans  of  the  flagged 
suspect  tissue  and  healthy  tissue  in  the  same  session. 
Access  to  data  from  6  patients  were  obtained  and  made 
available  for  analysis  in  this  study.  We  proposed  to  use  a 
combination  of  pattern  recognition  techniques,  including 
Artificial  Neural  Networks  (ANN) ,  to  develop  in  vivo 
methods  that  use  breast  P-31  MR  scans  (suspicious  and  non- 
suspicious  regions)  to  diagnose  potential  malignant 
tissue.  The  MR  scan  data  will  be  paired  with  the  known 
biopsy  results  to  create  a  supervised  training  set. 
Unfortunately  two  events  occurred  to  prevent  us  from 
completing  this  study[l-3]. 

Eirst,  originally  we  were  to  obtain  data  on  25 
patients.  However  before  all  25  patients  were  scanned, 
the  MRI  instrument  used  in  this  study  was  decommissioned 
for  replacement  by  a  more  powerful  and  updated  MRI .  We 
successfully  completed  MRI  P-31  scans  on  only  six 
patients.  The  transition  from  old  to  new  MRI  required 
approximately  4  months.  Second,  the  new  MRI  needed  to 
have  a  new  P-31  MRI  probe  constructed.  Such  a  task 
requires  expertise  in  probe  design.  The  expert  in  this 
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project  with  this  expertise.  Dr.  Diana  Lindquist,  took 
another  job  in  Cincinnati,  Ohio  before  she  was  able  to 
develop  a  new  probe  for  the  new  MRI . 

For  the  reasons  mentioned  previously,  it  is  highly 
desirable  to  develop  a  non-invasive  spectroscopic  and 
pattern  recognition  technique  that  can  detect  and  reliably 
interpret  images  or  spectral  data  from  small  volumes  of  the 
breast.  Due  to  the  pervasive  nature  of  breast  cancer  in 
society  today,  and  the  consequent  need  of  a  highly  accurate, 
early  diagnostic  tool,  this  is  a  very  timely  proposal  that 
could  have  a  significant  impact  on  women's  health.  The 
results  of  the  brain  scan  study  suggest  that  breast  tumor 
detection  may  be  possible  at  a  very  early  stage  with  this 
technique  given  sufficient  examples  in  the  training  sets 
[4]  .  There  is  a  high  degree  of  correlation  between  the 
chance  of  complete  eradication,  recovery,  and  survival  from 
cancer  and  how  early  it's  detected.  This  is  true  for  breast 
cancer  as  it  is  for  all  other  types  of  cancer. 

II  Methods : 


A.  Subjects: 

Seven  patients  were  chosen  because  they  had  a 
suspicious  mass  found  either  by  mammography  or  palpation 
and  were  scheduled  for  either  surgery  or  biopsy.  They  were 
excluded  if  they  had  had  a  previous  core  biopsy  from  the 
site.  The  local  institutional  review  board  approved  this 
study.  A  typical  MRI  instrument  containing  patient  is 
shown  in  Figure  la  with  corresponding  schematic  of 
important  components  shown  in  Figure  lb. 


SHIELD! 


MAGNET 


PATIENT 

TABLE 


RF  COIL 


GRADIENT  COILS 


MAGNET 


FILM 


GRADIENT 

AMP 


GRADIENT 
PULSE  PROG 


COMPUTER 


R 


RF 

DETECTOR 


DIGITIZER 


PULSE  PROG 


RF  SOURSE 


- RFAMP 


Figure  la:  Typical  MRI  Figure  lb:  Schematic  of 

instrument  with  a  patient  typical  MRI  instrument  used  . 

shown  being  readied  for  in  this  project. 
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MRI  scan. 

B.  MR  Data  Acquisition 

All  data  were  acquired  on  a  1.5  T  General  Electric 
Signa  Horizon  system  (General  Electric,  Milwaukee,  WI,  USA) 
using  a  home-built  solenoid  phosphorus  coil  placed  inside  a 
proton  transmit /receive  bilateral  breast  coil  (MRI  Devices, 
Milwaukee,  WI,  USA)  .  For  these  studies,  only  one  breast 
was  examined,  so  the  proton  coil  was  operated  as  either  a 
right-  or  left-only  coil. 

C.  Phosphorus -31  coil: 

The  phosphorus  coil  consisted  of  k-inch  wide  copper  tape 
wrapped  on  an  acrylic  tube  13.6  cm  in  diameter.  The  coil 
was  9  cm  tall  from  top  to  bottom.  The  acrylic  tube  was 
fitted  with  a  plastic  flange  to  allow  it  to  hang  from  the 
top  of  the  proton  breast  coil.  Variable  capacitors  were 
used  for  tuning  and  matching  the  coil,  which  was  done  for 
most  of  the  subjects  recruited.  A  balun  was  inserted  in 
the  coil  cable  to  minimize  effects  from  cable  currents. 
Figure  2  below  shows  a  typical  MRI  spectra  from  normal 
human  tissue  [1] . 


y-NTP 


ppm 


Figure  2.  Example  of  a  magnetic  resonance  spectrum  from 
normal  human  tissue. 

In  this  magnetic  resonance  spectrum  from  normal  human 
tissue  we  see  that  the  peak  area  is  proportional  to  amount 
of  metabolite.  Peaks  labeled  on  scan:  PME, 
phosphomonoesters  (mainly  phospholipid  precursors, 
principally  phosphocompounds  of  choline  and  ethanolamine. 


6 


with  some  contribution  from  sugar  phosphates)  ;  Pi, 
inorganic  phosphate  (product  of  ATP  hydrolysis) ;  PDE, 
phosphodiesters  (phospholipid  catabolites,  mainly 

glycerophosphocompounds  of  choline  and  ethanolamine,  with 
some  contribution  from  cell  membranes)  ;  and  y-,  a-,  (3- 

phosphates  of  nucleotide  triphosphates,  NTP  (high-energy 
phosphate  compounds)  .  By  convention,  the  [p-P]  NTP  peak  is 
taken  to  represent  ATP.  Phosphocreatine  peak  (not  labeled) 
is  at  zero  parts  per  million  (assumed  from  muscle 
contamination) .  Data  are  conventionally  presented  as  ratios 
of  related  peak  areas,  comprising  phosphoester  metabolites 
(PME/PDE)  and  energy  status  (ATP/Pi)  respectively. 
Alternatively,  individual  peak  areas  may  be  expressed  as  a 
function  of  total  visible  phosphate  (TP)  .  These  measures 
are  independent  of  the  volume  of  liver  from  which  the 
signals  are  obtained.  Liver  intracellular  pH  (pHic)  was 
derived  from  the  chemical  shift  difference  between  [a-P] 
NTP,  assigned  to  -7.5  ppm,  and  Pi,  according  to  the 
supplied  calibration  data,  ppm,  resonance  chemical  shift  in 
parts  per  million. 

The  phosphorus  coil  (see  for  example  Figures  3  and 
4)  [2]  was  centered  in  the  appropriate  side  of  the  proton 
coil  and  the  subject  was  positioned  prone  on  top  of  the 
coils  with  the  breast  to  be  examined  hanging  freely.  The 
contralateral  breast  was  compressed  against  the  chest  to 
minimize  artifacts  in  the  clinical  images.  The  patient  was 
then  centered  in  the  magnet  (see  Figure  la)  .  Pre-  and 
post-contrast  clinical  RoDEO  images  were  acquired  to 
localize  the  lesion  (parameters:  flip  angle:  45  degrees; 
repetition  time  20ms;  echo  time  5ms;  128x256  matrix;  slice 
thickness  1.4-1. 8  mm;  128  slices)  .  A  slice  containing  the 
lesion  was  chosen  and  used  as  the  center  point  for  the 
spectroscopy  exam.  The  spectroscopy  exam  used  a  phase- 
encoded  pulse  acquire  technique  (parameters:  25  degree  flip 
angle,  repetition  time  300ms,  6x6  matrix,  512  complex 

points,  2000  Hz  spectral  width,  75  NEX,  200  mm  FOV,  and  30 
mm  slice  thickness) .  Each  spectral  data  set  required  about 
30  minutes  to  acquire.  To  assess  the  reproducibility  of  the 
technique,  spectra  were  acquired  on  three  separate  days 
from  one  volunteer.  For  all  other  subjects,  spectra  were 
acquired  once. 
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^  High  Q 


Figure  3.  A  saddle-shaped  radio 
frequency  (RF)  coil  inside  a 
solenoid  magnet.  Bo=the  direction 
of  the  magnetic  field  of  the 
solenoidal  magnet,  Bi=the  direction 
of  the  magnetic  field  produced  by 
current  flow,  i,  in  the  RF  coil. 

Data  were  transferred  to  an  off-line  workstation  (SGI 
O2  system)  for  processing  and  analysis.  The  data  were 
processed  using  SA/GE  (Spectral  analysis  software.  General 
Electric,  Milwaukee,  WI,  USA)  .  All  data  were  given  10  Hz 
of  exponential  line-broadening  prior  to  Eourier 
transformation.  The  data  were  phased  automatically,  with 
some  manual  phasing  of  individual  spectra  if  needed.  The 
phased  spectra  were  overlaid  on  the  central  image  from  the 
clinical  exam  for  display  [5-7]  . 

Ill  Key  Research  Accomplishments : 

•  Ran  MRIs  on  six  patients  using  RoDEO  enhancements  to 
pinpoint  potential  breast  tissue  anomalies. 

•  Ran  36  voxels  of  data  for  P-31  spectra  on  all  6 
patients . 

•  Thirty-six,  180  cm^  voxels  were  obtained  and  the 
corresponding  P-31  spectra  were  obtained  for  each  MRI 
scan  of  each  patient. 

•  The  36  spectra  per  patient  were  examined  to  locate  the 

peaks  from  phosphomonoesters ,  phosphodiesters , 

phosphocreatine,  inorganic  phosphate,  and  the  ATP 
peaks . 

•  Analysis  was  performed  with  AI  software.  The  signal- 
to-noise  ratio  was  too  low  to  allow  the  unambiguous 
determination  of  either  the  location  or  intensity  of 
any  of  these  peaks  [6,7]  . 


Frequency 

Figure  4. Dependence  of 
bandwidth,  Af,  of  coil 
on  value  of  Q. 
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IV.  Reportable  Outcomes; 

The  following  are  examples  of  both  the  MRI  images  of  a 
breast  (Figure  5a  and  Figure  6a)  and  the  corresponding  P-31 
spectra  (Figure  5b  and  Figure  6b) of  the  lesion  which  was 
targeted.  The  "aqua"  spectra  seen  on  the  breast  image  is 
located  in  the  area  of  the  voxel  where  the  P-31  spectra  was 
obtained.  Spectra  were  obtained  for  both  normal  (benign) 
tissue  and  potentially  malignant  tissue. 


Figure  5  a;  MRI  image  of  breast  for 
Patient  1  (normal  tissue). 


Figure  5b;  P-31  speetra  of  breast  on  left 
for  Patient  1 . 


The  following  figures  (Figure  6a  and  6b)  show  the  MRI  image  and  eorresponding  P-31 
speetra  for  a  patient  with  a  diagnosed  malignant  tumor. 


Figure  6a;  MRI  image  of  breast  for 
Patient  2  (suspeeted 
malignant  area) 


Figure  6b;  P-31  speetra  of  breast  on  left 
for  Patient  2. 
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The  raw  data  from  the  phosphorus  scan,  as  previously 
stated,  was  transferred  to  an  SGI  O2  workstation  for  offline 
processing  using  the  GE  Spectral  Analysis  (SA/GE)  package. 
The  data  was  processed  to  achieve  the  best  resolution  and 
signal  to  noise  ratios,  which  required  zero-filling  and 
apodizing  the  raw  data  prior  to  Eourier  transformation. 
The  precise  processing  parameters  will  be  determined  from 
the  first  data  set  and  will  be  used  consistently  throughout 
this  project.  Initial  processing  parameters  will  be  zero¬ 
filling  once  with  10  Hz  exponential  line  broadening 
applied . 

V.  Conclusion  and  Siommary; 

With  breast  cancer  being  the  number  one  cancer  in  woman 
with  over  200,000  diagnosed  this  year,  and  more  than  1.6 
million  women  a  survivor  of  this  cancer,  we  strongly 
believe  that  the  procedure  outlined  in  this  research 
proposal  has  the  potential  to  provide  a  much  earlier 
diagnosis  which  should  reduce  the  need  of  biopsies,  and, 
potentially  provide  a  much  higher  cure  rate  due  to  this 
earlier  detection  and  diagnosis  [11,12]. 

Due  to  the  limited  number  of  patients  we  were  able  to 
recruit  for  this  project  it  was  not  possible  to  perform 
sufficient  artificial  neural  network  analysis.  We 
originally  wanted  to  study  at  least  25  patients,  of  which 
at  least  half,  or  approximately  10  to  12  were  diagnosed 
with  a  malignant  tumor  verified  using  a  biopsy.  However,  we 
were  only  able  to  recruit  6  patients  with  only  one  patient 
(patient  2)  having  a  confirmed  malignant  breast  tumor.  In 
addition,  the  Co-PI,  Dr.  Diana  Lindquist,  is  no  longer 
available  for  this  project.  She  has  accepted  as  of 
May, 2006  another  position  in  Cleveland,  Ohio.  Therefore  we 
will  not  be  able  to  perform  any  additional  MRI  scans  with 
the  accompanying  P-31  spectra  for  analysis. 

However,  in  discussions  with  Dr.  Lindquist  and  other 
Co-PIs  on  the  project,  we  will  begin  looking  at  in  vitro 
tissue  samples  which  have  undergone  a  biopsy  for  the 
presence  of  malignant  tissue.  We  will  perform  NMR  P-31 
spectral  analysis  on  these  tissue  samples  for  analysis  with 
ANN.  We  should  be  able  to  perform  analysis  on  several 
hundred  patient  samples  and  compare  these  malignant  samples 
with  non-malignant  samples. 
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VII  Appendix: 


Table  1 . 

spectral 
patient . 


The  following  data  is 
data  obtained  from  one 
This  is  only  one  data 


an  example  of  MRI  P-31 
run  of  one  voxel  from  one 
set  of  36  per  patient. 


-48.3372 

-910.782 

-3245.02 

-48.3372 

-910.782 

-3245.02 

-48.3372 

-910.782 

-3245.02 

-48.2428 
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-988.529 

-44.7497 

-755.977 

-988.529 
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-555.958 
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-935.479 

172.0492 

-42.0118 

-935.479 

172.0492 
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-711.475 

1285.723 
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-711.475 

1285.723 

-41.9174 

-711.475 
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1887.929 

-41.7286 

1118.031 
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-40.8789 
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436.7585 

800.9824 

-40.7845 

436.7585 

800.9824 

-40.6901 

655.709 

1258.421 

-40.6901 

655.709 

1258.421 

-40.6901 

655.709 

1258.421 

-40.5957 

1265.602 

1291.614 

-40.5957 

1265.602 

1291.614 

-40.5957 

1265.602 

1291.614 

-40.5013 

1546.083 

990.8132 

-40.5013 

1546.083 

990.8132 

-40.5013 

1546.083 

990.8132 

-40.4069 

1661.777 

774.5599 

-40.4069 

1661.777 

774.5599 

-40.4069 

1661.777 

774.5599 

-40.3125 

1750.604 

542.1177 

-40.3125 

1750.604 

542.1177 

-40.3125 

1750.604 

542.1177 

-40.2181 

1698.236 

223.1239 

-40.2181 

1698.236 

223.1239 

-40.2181 

1698.236 

223.1239 

-40.1237 

1400.869 

47.88004 

-40.1237 

1400.869 

47.88004 

-40.1237 

1400.869 

47.88004 

-40.0292 

1095.401 

106.7472 

-40.0292 

1095.401 

106.7472 

-40.0292 

1095.401 

106.7472 

-39.9348 

785.3507 

344.7289 

-39.9348 

785.3507 

344.7289 

-39.9348 

785.3507 

344.7289 

-39.8404 

581.1136 

966.6189 

-39.8404 

581.1136 

966.6189 

-39.8404 

581.1136 

966.6189 

-39.746 

1094.592 

1895.979 

-39.746 

1094.592 

1895.979 

-39.746 

1094.592 

1895.979 

-39.6516 

2400.454 

2055.838 

-39.6516 

2400.454 

2055.838 

-39.6516 

2400.454 

2055.838 

-39.5572 

3264.755 

1286.623 

-39.5572 

3264.755 

1286.623 

-39.5572 

3264.755 

1286.623 

-39.4628 

3576.869 

486.9435 

-39.4628 

3576.869 

486.9435 

-39.4628 

3576.869 

486.9435 

-39.3684 

3654.662 

-417.049 

-39.3684 

3654.662 

-417.049 

-39.3684 

3654.662 

-417.049 

-39.274 

3095.014 

-1304.43 

-39.274 

3095.014 

-1304.43 

-39.274 

3095.014 

-1304.43 

-39.1796 

2165.288 

-1392.04 

-39.1796 

2165.288 

-1392.04 

-39.1796 

2165.288 

-1392.04 
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-39.0852 

1868.401 

-851.317 

-39.0852 

1868.401 

-851.317 

-39.0852 

1868.401 

-851.317 

-38.9907 

2328.221 

-664.094 

-38.9907 

2328.221 

-664.094 

-38.9907 

2328.221 

-664.094 

-38.8963 

2662.348 

-1368.09 

-38.8963 

2662.348 

-1368.09 

-38.8963 

2662.348 

-1368.09 

-38.8019 

2035.069 

-2134.88 

-38.8019 

2035.069 

-2134.88 

-38.8019 

2035.069 

-2134.88 

-38.7075 

1126.071 

-2161.85 

-38.7075 

1126.071 

-2161.85 

-38.7075 

1126.071 

-2161.85 

-38.6131 

427.0704 

-1735.19 

-38.6131 

427.0704 

-1735.19 

-38.6131 

427.0704 

-1735.19 

-38.5187 

117.3701 

-998.099 

-38.5187 

117.3701 

-998.099 

-38.5187 

117.3701 

-998.099 

-38.4243 

270.3219 

-340.054 

-38.4243 

270.3219 

-340.054 

-38.4243 

270.3219 

-340.054 

-38.3299 

605.017 

152.7145 

-38.3299 

605.017 

152.7145 

-38.3299 

605.017 

152.7145 

-38.2355 

1298.704 

543.0816 

-38.2355 

1298.704 

543.0816 

-38.2355 

1298.704 

543.0816 

-38.1411 

2218.775 

161.337 

-38.1411 

2218.775 

161.337 

-38.1411 

2218.775 

161.337 

-38.0467 

2298.62 

-648.663 

-38.0467 

2298.62 

-648.663 

-38.0467 

2298.62 

-648.663 

-37.9523 

2065.599 

-883.108 

-37.9523 

2065.599 

-883.108 

-37.9523 

2065.599 

-883.108 

-37.8578 

2018.563 

-1031.5 

-37.8578 

2018.563 

-1031.5 

-37.8578 

2018.563 

-1031.5 

-37.7634 

1994.053 

-1019.03 

-37.7634 

1994.053 

-1019.03 

-37.7634 

1994.053 

-1019.03 

-37.669 

2434.039 

-1078.69 

-37.669 

2434.039 

-1078.69 

-37.669 

2434.039 

-1078.69 

-37.5746 

2926.522 

-1836.27 

-37.5746 

2926.522 

-1836.27 

-37.5746 

2926.522 

-1836.27 

-37.4802 

2797.235 

-3020.37 

-37.4802 

2797.235 

-3020.37 

-37.4802 

2797.235 

-3020.37 

-37.3858 

1955.65 

-4154.11 

-37.3858 

1955.65 

-4154.11 

-37.3858 

1955.65 

-4154.11 

-37.2914 

404.4081 

-4773.69 

-37.2914 

404.4081 

-4773.69 

-37.2914 

404.4081 

-4773.69 

-37.197 

-1253.91 

-4391.2 

-37.197 

-1253.91 

-4391.2 

-37.197 

-1253.91 

-4391.2 

-37.1026 

-2350.22 

-3401.23 

-37.1026 

-2350.22 

-3401.23 

-37.1026 

-2350.22 

-3401.23 

-34.7424 

3228.471 

-2040.56 

-34.7424 

3228.471 

-2040.56 

-34.7424 

3228.471 

-2040.56 

-34.648 

3100.886 

-3279.93 

-34.648 

3100.886 

-3279.93 

-34.648 

3100.886 

-3279.93 

-34.5535 

1600.654 

-4325.02 

-34.5535 

1600.654 

-4325.02 

-34.5535 

1600.654 

-4325.02 

-34.4591 

-486.409 

-3841.55 

-34.4591 

-486.409 

-3841.55 

-34.4591 

-486.409 

-3841.55 

-34.3647 

-1738.76 

-1897.63 

-34.3647 

-1738.76 

-1897.63 

-34.3647 

-1738.76 

-1897.63 

-34.2703 

-1260.44 

472.7027 

-34.2703 

-1260.44 

472.7027 

-34.2703 

-1260.44 

472.7027 

-34.1759 

703.5732 

1766.938 

-34.1759 

703.5732 

1766.938 

-34.1759 

703.5732 

1766.938 

-34.0815 

2688.591 

1443.014 

-34.0815 

2688.591 

1443.014 

-34.0815 

2688.591 

1443.014 

-33.9871 

3933.663 

426.366 

-33.9871 

3933.663 

426.366 

-33.9871 

3933.663 

426.366 

-33.8927 

4650.686 

-1056.43 

-33.8927 

4650.686 

-1056.43 

-33.8927 

4650.686 

-1056.43 

-33.7983 

4258.882 

-2612.09 

-33.7983 

4258.882 

-2612.09 

-33.7983 

4258.882 

-2612.09 

-33.7039 

3496.437 

-3301.98 

-33.7039 

3496.437 

-3301.98 

-33.7039 

3496.437 

-3301.98 

-33.6095 

2949.12 

-3946.55 

-33.6095 

2949.12 

-3946.55 

-33.6095 

2949.12 

-3946.55 

-33.5151 

1877.297 

-4492.99 

-33.5151 

1877.297 

-4492.99 

-33.5151 

1877.297 

-4492.99 

-33.4206 

579.0143 

-4356.7 

-33.4206 

579.0143 

-4356.7 

-33.4206 

579.0143 

-4356.7 

-33.3262 

-573.111 

-3364.94 

-33.3262 

-573.111 

-3364.94 

-33.3262 

-573.111 

-3364.94 

-33.2318 

-498.53 

-1619.73 

-33.2318 

-498.53 

-1619.73 

-33.2318 

-498.53 

-1619.73 

-33.1374 

984.8907 

-942.968 

-33.1374 

984.8907 

-942.968 

-33.1374 

984.8907 

-942.968 

-33.043 

1890.062 

-1726.66 

-33.043 

1890.062 

-1726.66 

-33.043 

1890.062 

-1726.66 

-32.9486 

1781.755 

-2332.74 

-32.9486 

1781.755 

-2332.74 

-32.9486 

1781.755 

-2332.74 

-32.8542 

1768.243 

-2494.29 

-32.8542 

1768.243 

-2494.29 

-32.8542 

1768.243 

-2494.29 

-32.7598 

1813.553 

-2827.02 

-32.7598 

1813.553 

-2827.02 

-32.7598 

1813.553 

-2827.02 

-32.6654 

1756.955 

-3052.58 

-32.6654 

1756.955 

-3052.58 

-32.6654 

1756.955 

-3052.58 
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